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Abstract 

An evaluation of the effective QED coupling at the scale Mz is presented. It em- 
ploys the predictions of perturbative QCD for the cross section of electron positron 
annihilation into hadrons up to order a^, including the full quark mass depen- 
dence, and of order in the high energy region. This allows to predict the input 
for the dispersion relations over a large part of the integration region. The per- 
turbative piece is combined with data for the lower energies and the heavy quark 
thresholds. The result for the hadronic contribution to the running of the coupling 
Aa[fjj(M|) = (277.5±1.7)xl0-^ leads to (q(M|))-i = 128.927±0.023. Compared 
to previous analyses the uncertainty is thus significantly reduced. 



The evolution of the electromagnetic coupling from its definition at vanishing mo- 
mentum transfer to its value at high energies constitutes the dominant part of radiative 
corrections for electroweak observables. The accurate determination of a(M|) is thus 
essential for any precise test of the theory. At the same time the indirect determina- 
tion of the masses of heavy, hitherto unobserved particles, e.g. the Higgs boson or SUSY 
particles, depends critically on this parameter. Of particular importance in this context 
is the hadronic vacuum polarization. It is nearly as large as the leptonic contribution 
but can only be related through dispersion relations to the cross section for hadron pro- 
duction in electron positron annihilation, or more conveniently to the familiar R ratio. 
The integrand can thus be obtained from data, phenomenological models and/or per- 
turbative QCD (pQCD), whenever apphcable. A detailed analysis based on data and 
employing pQCD above 40 GeV has been performed in [|^ and their result has been 
confirmed by subsequent studies [Q, Q following very similar strategies. Alternatively 
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one may employ pQCD also for lower energies, eventually as low as 2 GeV as long as 
one stays away from the quark thresholds. A first step in this direction has been made 
in 0. There, however, only the massless approximation was employed for the normal- 
ization of the data. Recently the 0{a'^) corrections for R, including the full quark mass 
dependence, became available ||, |^, |^ which allows to extend pQCD down relatively close 
to the respective thresholds for charm and bottom production. These results have been 
used in P] to evaluate the perturbative contribution to the vacuum polarization. In this 
work we complete the evaluation by incorporating those contributions which cannot be 
obtained from pQCD: from the low energy region below roughly 2 GeV, the charmonium 
and bottomonium resonances and from the charm threshold. The (re)normalization of 
the data from the PLUTO 0, DASP lTl|, and MARK I collaborations ^ will be based 
on the requirement that they agree with pQCD for y/s < 3.7 GeV and y/s > 5 GeV. 



A similar analysis has been recently performed in |jT3[ which is also based on pQCD, in 
particular the results of [0 and which, furthermore, provides additional justification for 
the applicability of pQCD at very low energies around 2 GeV. We will comment on the 
differences between this analysis and ||TB[, whenever appropriate, below. 

Let us briefly describe the theoretical input for our evaluation. The hadronic contri- 
bution Aa[|^j to the running from the static limit to Mz is given by 



Svr J4m2 sis — My — It) 



where the superscript "(5)" indicates that the top quark is not included in the integral. 
For the analysis we use a = a(0) = 1/137.0359895 and Mz = 91.187 GeV. It leads, after 
resummation of the leading logarithms, to the following shift 



a(0) 

1 - Aaicp(s) - Aajf^dl-s) - Aatop(s) 



= , . . — (2) 



The quantity R{s) can be experimentally determined through a measurement of the total 
cross section for electron positron annihilation into hadrons. From the theoretical side it 
is defined through the absorptive part of the electromagnetic current correlator 

(-gV + ^M^O = ^ldxe'^^{0\Tj^ix)UO)\0), (3) 

R{s) = 12n lmn{q^ = s + ie) . (4) 

It can be calculated in the framework of pQCD up to order if quark masses are 
neglected |T^ and up to 0{al) with full quark mass dependence 0, |[. In this work 
pQCD will be assumed to be valid above 1.8 GeV (alternatively 2.125 GeV). In view of 
the validity of pQCD in tau lepton decays not only for the total rate but also for the 
spectral function toward the upper end [|15] the substitution of inprecise data by pQCD 
seems justified. Additional support for this approach can also be drawn from the analysis 



of data for R{s) below 1.8 GeV |T3|. pQCD has also been used in [|T6| in the present 



context even down to 1.4 GeV. The specific choice of 1.8 GeV (or 2.125 GeV) for the 
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matching between data and theory is dictated by the available data analysis ^ which 
we adopt for the present purpose. 

Also excluded from this theory driven treatment are the threshold region for charmed 
mesons (the interval from 3.7 GeV to 5 GeV) and, similarly, for bottom mesons (10.5 GeV 
to 11.2 GeV), and the narrow charmonium and bottomonium resonances, where we use 
the currently available data. Perturbative QCD is even applicable in the charm and 
bottom threshold regions, as far as the light quark contributions are concerned. In the 
bottom threshold region we will therefore use data for the bottom contribution only. 

In the perturbative regions one receives contributions from light {u, d and s) quarks 
whose masses are neglected throughout, and from massive quarks which demand a more 
refined treatment. Below the charm threshold the light quark contributions are evaluated 
in order plus terms of order a^s/(4M^) from virtual massive quark loops. Above 
5 GeV the full Mc dependence is taken into account up to order a^, and in addition the 
dominant cubic terms in the strong coupling are incorporated, as well as the corrections 
from virtual bottom quark loops of order aj. s/ {4M^). Above 11.2 GeV the same formalism 
is applied to the massive bottom quarks and charmed mass effects are taken into account 
through their leading contributions in an M^/s expansion. All formulae are available for 
arbitrary renormalization scale fi which allows to test the scale dependence of the final 
answer. This will be used to estimate the theoretical uncertainties from uncalculated 
higher orders. Matching of between the treatment with n/ = 3, 4 and 5 flavours is 
performed at the respective threshold values. The influence of the small 0{al) singlet 
piece which prevents a clear separation of contributions from different quark species can 
safely be ignored for the present purpose. The details of the formalism can be found in 

In Tab. |I|, adopted from 0, the perturbative hadronic contributions are hsted sepa- 
rately for a variety of energy intervals. As our default values we adopt /i^ = s, a^^)(M|) = 
0.118, Mc = 1.6 GeV and = 4.7 GeV. In separate columns we list the variations with 
a change in the renormalization scale, the strong coupling constant and the quark masses: 

5as = ±0.003, 5Mc = ±0.2 GeV, 5Mb = ±0.3 GeV. (5) 

In principle the theoretical tools are available to include in the perturbative analysis 
the QED corrections of order a and even aag- The relative size of the dominant correction 
is estimated as 



(0.6-0.7) X 10"% (6) 



and is also included in Tab. [|. 

Perturbative QCD is clearly inapplicable in the charm threshold region between 3.7 and 
5 GeV where rapid variations of the cross section are observed. Data have been taken more 



than 15 years ago by the PLUTO |g, DASP [y], and MARK I collaborations [|T2[. The 
systematic errors of 10 to 20 % exceed the statistical errors significantly and are reflected 
in a sizeable spread of the experimental results. To arrive at a reliable evaluation of the 
charm contribution from this region, we adjust the normalization of the data (for each 
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Energy range (GeV) 


central value 


6n 


6as 


6Mc 


6Mb 


1.800 - 2.125 


5.67 


0.22 


0.04 


0.00 


0.00 


2.125 - 3.000 


11.66 


0.21 


0.06 


0.01 


0.00 


3.000 - 3.700 


7.03 


0.06 


0.03 


0.00 


0.00 


1.800 - 3.700 


24.36 


0.48 


0.13 


0.01 


0.01 


5.000 - 5.500 


5.44 


0.03 


0.03 


0.06 


0.00 


5.500 - 6.000 


4.93 


0.03 


0.02 


0.04 


0.00 


6.000 - 9.460 


25.45 


0.11 


0.08 


0.10 


0.00 


9.460 - 10.520 


5.90 


0.02 


0.01 


0.01 


0.00 


10.520 - 11.200 


3.48 


0.01 


0.01 


0.00 


0.00 


5.000 - 11.200 


45.20 


0.19 


0.15 


0.21 


0.01 


(without bb) 












11.200 - 11.500 


1.63 


0.00 


0.01 


0.00 


0.00 


11.500 - 12.000 


2.62 


0.00 


0.01 


0.00 


0.00 


12.000 - 13.000 


4.93 


0.01 


0.01 


0.00 


0.00 


13.000 - 40.000 


72.92 


0.08 


0.12 


0.02 


0.02 


12.000 - 40.000 


77.85 


0.09 


0.14 


0.02 


0.02 


40.000 - oo 


42.67 


0.03 


0.06 


0.00 


0.00 


11.200 - oo 


124.77 


0.12 


0.21 


0.03 


0.02 


1.8 - oo (pQCD) 


194.33 


0.79 


0.49 


0.24 


0.03 


QED 


0.11 











Table 1: Contributions to Aq;[^J^(M|) (in units of 10""^) from the energy- 
regions where pQCD is used (adopted from 0). For the QED corrections the 
same intervals have been chosen. For the variation of as(M|), Mc and Mb 
Eqs. (^ have been used, fi has been varied between y/s/2 and 2y/s. 



experiment individually) to the theoretical predictions at the upper and lower endpoints 
as follows: Data below and up to 3.7 GeV are combined to determine the factor n_ which 
characterizes the mismatch between theory and experiment below threshold 




\RpQCd{s) 



and an averaged experimental R value just below threshold 

i?_ = n_/?pQCD((3.7GeV)2). 
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PLUTO 


DASP 


MARKl 


Interval below (GeV) 
ri- 

R- 


3.6000 - 3.6600 
1.04 ±0.01 
2.25 ±0.03 


3.6025 - 3.6500 
1.06 ±0.02 
2.29 ±0.04 


3.0000 - 3.6500 
1.18 ±0.04 
2.55 ±0.08 


Interval above (GeV) 
n+ 

R+ 


4.9800 - 4.9800 
1.04 ±0.01 
3.85 ± 0.04 


5.0000 - 5.1950 
1.15 ±0.01 
4.27 ±0.04 


5.1000 - 6.0000 
1.12 ±0.01 
4.14 ±0.05 


Aaif (M|) X 10^ (Model 1) 
Aaif (M|) X 10^ (Model 2) 


15.65 ±0.19 
15.64 ±0.16 


15.26 ±0.25 
15.68 ±0.30 


16.22 ±0.32 
15.83 ±0.36 



Table 2: Contribution to A«[^j^{j(M|) from the energy interval 3.7 to 5.0 GeV. 



In a similar way r2+ and R+ are derived from the data around and above 5 GeV. The 
normalization factors and the combined R values as given in Tab. ^ are consistent with 
the systematic errors quoted by the experiments. To account for the difference between 
n_ and rtj^ even within one experiment, two models are used for the interpolation into the 
interior of the interval. Model 1: The difference is due to the different efficiencies for final 
states with and without charmed mesons. Model 2: The difference is due to a linear s 
dependence of the experimental normalization and thus reflected in a linear s dependence 
of the quantity l/n{s) in the threshold region. 

The average of the two slightly different results is taken as central value and the three 
experiments are then assumed to be uncorrelated. The typical spread of ±0.2 is taken as 
systematical uncertainty which is combined linearly with the statistical error for which 
we take the maximum of Model 1 and Model 2. The combined result is thus given by 

Aaf^Ml) = (15.67 ± 0.34) X 10'^ (9) 

A similar approach has been adopted in [Q. In 0, however, only MARK I data were 
employed (with DASP, PLUTO and Crystal Ball data used for cross checks), an energy 
independent correction factor was chosen, and the pQCD prediction for massless quarks 
was used for the calibration below charm threshold. 

For the three lowest and narrow charmonium resonances we use the narrow width 
approximation: 

with (a/a(M|))^ = 0.96. Mr is the mass of the resonance and Fee the partial width 
into electrons. The errors from the three charmonium resonances are added linearly. The 



result, given in Tab. I, differs by about 0.7 x 10~* from |I| 
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Input 


energy region 


Aa(5) X 10^ 


low energy data IS] 


2m^ - 1.8 GeV 


56.90 ± 1.10 


narrow charmonium resonances 


J/^,^'(25),'f(3770) 


9.24 ± 0.74 


"normalized" data 


3.7 - 5.0 GeV 


15.67 ±0.34 


T resonances 


T{1S) - T(11.019) 


1.17 ±0.09 


interpolation of bb 


11.075 - 11.2 


0.03 ± 0.03 


pQCD (and QED) 


1.8 - oo 


194.45 ± 0.96 




total 


277.45 ± 1.68 



Table 3: Contributions to Aal^J^^{M'^) from different energy regions. 



The contributions from the three lowest T resonances are evaluated through Eq. ([T0|) 
with {a/a{My)y = 0.93. For the bottom threshold an approach similar to the one for 
charm could be employed. However, the bb cross section between 10.5 and 11.075 GeV is 
saturated by the three T resonances at 10.580 GeV, 10.865 GeV and 11.019 GeV. (The 
result for the six T resonances as listed in Tab. |^ differs from |13[ by 0.2 x 10~^.) For 



energies above 11.2 GeV the perturbative treatment seems adequate. Between 11.075 and 
11.2 GeV a linear increase from zero to the perturbative value is assumed, and the error 
is conservatively taken to be equal to this value. 

For the low energy region up to 1.8 GeV we use the value (56.90±1.10) x 10~^ |T^. The 



individual results for the different regions are listed in Tab. ^. Combining the experimental 
errors, those from q;^^^(M|), the quark masses and the theoretical error in quadrature, we 
find 

Aaa(M|) = (277.45 ± 1.68) X 10-^ (11) 

as our main result. Alternatively we could have used the more conservative analysis of 0] 
for the region up to 2.125 GeV with a contribution of (63.42 ± 2.59) x lO"'' and pQCD 
only above 2.125 GeV. The resuh of this approach, Aa[,5d(M|) = (278.30 ± 2.82) x 10~\ 
would differ slightly in the central value and significantly in the size of the error. 

Frequently the contribution from the top quark is added to this value. Using the 
three-loop QCD corrections 



'77-T72 S 1 + 5.062^^^ ± 28.220 ± 9.702 In 

45 IT TT V 




one obtains 



0.1071 ± 0.8315 ' ' + 6.924 ± 1.594 In ^ 
TT V Mi 



Aatop(M|) = (-0.70 ± 0.05) X 10-\ (13) 
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where we have used Mt = 175.6 ±5.5 GeV. For convenience Eq. ( p!2[ ) is expressed in terms 
of a^^^fi^). For the numerical evaluation a^^''(M|) = 0.118 has been chosen. 

A major contribution to the vacuum polarization originates from the leptons. The 
dominant term is given by 



/ , ^9 N " / 5 1 My mj ( mf 




z_ 

2 

2n , /n f^.3\ 



+ Aaiep,2i(Mi) ± O [a' 
314.19 X 10"^ ± Aaiep,2i(M|) ± O (a^) . (14) 



The two-loop correction 

A„„,,,(M|) = (-) E f__ + C(3) + -ln-| + 3^ln-f + 0(^jj 

(15) 

leads to a shift which could in principle become relevant in forthcoming precision studies. 
For the combined result we thus obtain 

(a(M|))"^ = 128.927 ±0.023, (16) 




if we use the more optimistic analysis |113| for the region below 1.8 GeV, and alternatively 
(a(M|))"^ = 128.916 ± 0.039 if we employ the analysis from § below 2.125 GeV. In 
Tab. I our result for Aajj^^^(M|) is compared to earlier evaluations. Our uncertainty is 
only a quarter of the one from the analysis of |jl| or Q based on data only — at the price 
of a more pronounced dependence on pQCD at relatively low energies. The reduction 
of the error by a factor 1.5 compared to |0 is to a large extend a consequence of our 
different treatment of the charm threshold. The shift of the central value by —1.0 x 10~^ 



compared to [13] is mainly due to different values for the charmonium and bottomonium 
contributions and our treatment of the charm threshold. In the prediction for a(M|) this 
is partly compensated by our inclusion of the leptonic two-loop contribution of 0.8 x 10~^. 

Summary: The effective fine structure constant at Mz has been evaluated with input 
from pQCD over most of the integration region. A detailed analysis of the theoretical 
uncertainties has been performed. The two-loop piece for leptons has been included. In 
comparison with earlier results based on the analysis of data a slight shift of the central 
value and a significant reduction of the error has been obtained. 
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Aa^l^l) X 


Reference 


273.2 ±4.2 
280 ± 7 
280 ± 7 
275.2 ±4.6 
281.7 ±6.2 

278.4 ± 2.6* 

277.5 ± 1.7 


Martin et al. '95 
[H, Eidelman et al. '95 
1], Burkhardt et al. '95 
1], Swartz '96 
1], Alemany et al. '97 
[13|, Davier et al. '97 
this work 



Table 4: Comparison of different evaluations of Aa|j2i(^I)- (*^cttop(Af|) 
subtracted.) 
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